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Extended instantons generated on the lattice 

E. MendeP* and G. Nolte ^ 

'^FB Physik, Carl von Ossietzky Universitat Oldenburg 
26111 Oldenburg, Germany 

We have been able to observe directly extended instantons on the lattice, with a new method that does not 
require dislocations to measure them, and where we do not perform cooling. We showed, based on the simple 
Abelian Higgs model in 1 + 1 dim., that one can extract the instanton and anti-instanton density and their size, 
by measuring the topological charge, Q„ , on sub- volumes v larger than the instanton sizes, but smaller than the 
periodic lattice of size V. We are working on the generalization for non-abelian models. 



1. Introduction 

Instantons ||l|,^ are expected to play an im- 
portant role in such diverse phenomena as net 
Baryon number generation, the U{1) problem or 
even quark confinement. In order to extract their 
Quantum contribution to the path integral it is 
desirable to use a lattice formulation. We will 
propose a method that enables the observation of 
these topological non trivial configurations, even 
for relatively smooth configurations that don't 
present lattice artifacts called dislocations. (The 
topological charge is usually only measured on 
the whole volume V where, as a consequence of 
periodic boundary conditions (PBC), we can only 
see a net charge if we have dislocations). Further- 
more, we will show that one can extract informa- 
tion on the instanton density pi and size Lj , with- 
out the need of cooling to suppress the quan- 
tum fluctuations. Note that one will still need 
improved actions Q or cooling, for 4 dimensional 
gauge theories, in order to avoid small instantons 
comparable with the lattice spacing a. If eventu- 
ally "perfect actions" are found that would avoid 
dislocations, our method could be very useful to 
enable the extraction of topological information. 

The method jsj, which we present here, con- 
sists in measuring the topological charge on sub- 
volumes V (of V), larger than or of the order of 
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the relevant instantons (I, AI). Let us first assume 
that one can reach a coupling regime where the 
gauge configurations are smooth enough so that 
there are no individual sites where the topologi- 
cal density gets even close to ±1/2. This can be 
reached for the abelian Higgs model in 1 -I- 1 di- 
mensions at weak coupling. Then one has safely 
configurations without dislocations and if we take 
PBC the total geometric topological charge is ex- 
actly zero. This just means that the charge for a 
number of I and AI has to add to zero. 

We will use as our main observable, the proba- 
bility distribution to get some topological charge 
Q in sub- volumes w, Pv{Q)- For very large v, this 
distribution will have lumps peaked at integer val- 
ues, representing a net number of I minus AI in 
V. As we reduce the volume v to sizes smaller 
than the instanton, we can only see a distribution 
peaked at Q = as we have just pieces of instan- 
tons in V or quantum fluctuations. Interestingly, 
these are Q fluctuations in regions of stable topo- 
logical vacua and therefore tend to cancel out over 
fairly short distances. Due to their topological 
character the total dispersion a, that they pro- 
duce in Q, only grows as with s being the free 
surface. For larger volumes, of the order or some- 
what larger than the instantons, we start getting 
relevant contributions to Pv{Q) from instantons, 
where the charge adds coherently to ±1 (if they 
are fully in v) in contrast to the fluctuations. 

The idea is then to look at scales v, where: 
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so that we have mainly suppressed the fluctua- 
tions, and can study the ensemble of I and AI 
with effectively almost free spatial BC. 

For the abelian Higgs model on fine lattices, 
we have found clearly identifiable extended in- 
stantons over several lattice spacings, while the 
charge over the whole V stays always zero. As- 
suming an almost free ensemble of I and AI, we 
have been able to extract an instanton average 
density pi and size Lj (also as function of temper- 
ature, T = 1/Nt). The density p/ which can also 
be extracted from the second moment of Pv{Q) 
for large v, falls steeply at higher T as the in- 
stantons get squeezed in the time direction. We 
have also studied ||^ the evolution of the posi- 
tion of the instantons in Monte Carlo time and 
found that they get created locally as I - AI pairs 
and annihilate later when they encounter another 
pairing possibility. 

As for some models it is hard to avoid disloca- 
tions, we have analyzed for our case what hap- 
pens if one takes rougher lattices where disloca- 
tions start to appear. Only the net number N of 
dislocations matters (pairs can be gauged away), 
effectively behaving as modified boundary condi- 
tions with N instantons. We have found in this 
coupling regime, that the probability to be in a 
sector with TV net instantons in the full volume 
gives a compatible distribution to the case with- 
out dislocations in sub-volumes of the same size. 

The previous discussion has been done using 
the geometrical definition for the topological den- 
sity 1^,0 but we have also investigated the good- 
ness of the field theoretic definition which is easier 
to implement ||^ and gives the right naive con- 
tinuum limit, but does have only approximately 
right addition properties and so does not give ex- 
actly integers on the full volume V. 

2. Topological Q for Abelian Higgs model 

The lattice formulation of the abelian Higgs 
model in 1-1-1 dim., which is one of the simplest 
allowing non-trivial topology, has the action Sn'- 

A($:$„-l)'-/3i?e([/„)-2Ki?e($;f/„,^$„+^)(l) 

Here $„ is the Higgs field, Un,^ = exp(z6'„.p) are 
the links and C/„ = HnJ7„^^ the 2-d plaquettes. 



In the continuum the topological charge in v is 
Qi, = T~ / '^t^iyFfj.f ^ ds nf^Ai_, . (2) 

so that on the whole V it is always zero for PBC. 

On the lattice the field-theoretic definition, giv- 
ing the right naive continuum limit, 

has the disadvantage of not being a proper topo- 
logical density, not satisfying exactly the second 
equality in eq.(^) and so it isn't strictly in V. 

The geometrical charge definition, which is the 
simplest example of the fibre bundle construction 
used in the non-abelian case is 

This may still violate eq.(||) by an integer 
amount in cases where the phases in the links of 
a plaquette add up to more than ±7r, this value 
being reduced by ±27r, by the log to its princi- 
pal branch. These lattice artifacts, called dislo- 
cations, are the reason for obtaining net charge on 
a periodic lattice. These artifacts can be avoided 
on fine enough lattices at least for this model. 

In order to keep track of these dislocations we 
will directly work with the phases 6n,p, of the links 
as the fundamental variables. For the plaquettes. 
On ^ 6n,x+0n+x,f-dn+f,x-6n,f- We then define 

--^ ^Y.\-On] , (5) 

where the brackets shift On by an integer multiple 
of 27r into the interval (— tt, tt]. This definition is 
clearly equivalent to (^, however, using 6'„ as the 
fundamental variable we are now able to locate 
dislocations, namely places for which |0„| > tt. 
Note that in our Metropolis algorithm we update 
the On,^i by small changes from the old values, oth- 
erwise we could trivially generate dislocations by 
adding 2-k to a link. We find that if we choose pa- 
rameters corresponding to a fine enough lattice, 
so that the physical instantons are much larger 
than a, we can avoid obtaining any dislocations 
at all over the whole Monte Carlo run. Most sim- 
ulations have been done then around the "scaling 
region " §, (A = .2, k ^ .37, /? > 7.). 
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Top. charge distr. for loops 6*160 in 6*330 latt. 
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Charge Q 

Figure 1. Probability to find a charge Q in a 
volume V, half of V. Note the data agreement 
with the model over 3 orders of magnitude. The 
rise at integers is due to instantons fully in v. 

3. Identifying Instantons 

In order to be able to identify instantons 
clearly, we have chosen lattices with spatially 
very elongated V and with subvolumes v which 
wrap around the r axis in order to suppress the 
charge fluctuations along that long "boundary" . 
The finite temperature also serves to get a di- 
luter system of instantons. We took for example 
F = 320 X 6 and w = (20 - 160) x 6 and could 
clearly see the creation of I - AI pairs. 

The charge probability distribution Pv{Q) is 
shown in fig. 1, with clear peaks at integer charge. 
For a quantitative analysis to extract pi, Lj and 
the dispersion in the fluctuations a, we developed 
a simple model of almost noninteracting instan- 
tons. Assuming a number density pi = pAi, the 
probability to have n pairs in the total V is 



PnipV) 



1 



i!2 N 



(6) 



where N is the normalization. This p„ has to be 
folded with the number of ways to get a charge Q 
in V, assuming that the n instantons of size Lj can 
be located anywhere in V. With this distribution 
fniQ) in a given sector, we get 



This ideal distribution can be smeared by a gaus- 
sian with a width a due to quantum fluctuations. 
The best fit for the parameters in lattice units is 
p w .0011, L/ « 28 and (7 « .27. The topological 
susceptibility can also be modeled |^ for large V 
being < >« 2vpi. We have checked compati- 
bility for various v, and calculated at some T. 

4. Beyond optimal conditions 

As for some models we cannot reach couplings 
where we can suppress dislocations, it is interest- 
ing to compare ioi a, (3 — 6.3, where they just 
start to appear, the PydiQ) for sectors with net 
number of dislocations (and net instantons) with 
the Pv{Q) without dislocations with v = Vd- The 
agreement is shown in Tabic 1 , telling us that the 
net instantons come with the right S and entropy. 
In this sense the traditional way to measure topol- 
ogy is fine, at least for large instantons. 

We also used the field theoretic and found 
good agreement |^ with for fine lattices. For 
SU{2) we have observed much smaller fluctua- 
tions for if, instead of symmetrizing loops as 
in Ref.^, we reorder terms (same Q on full V) 
taking for each hypercube all loops wrapping it. 

Table 1 

Relative weights of sectors with charge |n|. 



n 





1 


2 


3 


with disloc. 


1. 


.51(4) 


.11(1.5) 


.005(2) 


no disloc. 


1. 


.53(2) 


.10(1) 


.007(3) 
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